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We investigate the edge reconstruction phenomenon beheved to occur in quantum dots in the 
quantum Hall regime when the filling fraction is v ^ 1. Our approach involves the examination of 
large dots (< 40 electrons) using a partial diagonalization technique in which the occupancies of the 
deep interior orbitals are frozen. To interpret the results of this calculation, we evaluate the overlap 
between the diagonalized ground state and a set of trial wavefunctions which we call projected 
necklace (PN) states. A PN state is simply the angular momentum projection of a maximum 
density droplet surrounded by a ring of localized electrons. Our calculations reveal that PN states 
have up to 99% overlap with the diagonalized ground states, and are lower in energy than the states 
identified in Chamon and Wen's study of the edge reconstruction. 

PACS numbers: 71.10.-w, 73.23.Hk, 73.40. Hm, 73.61.-r 



I. INTRODUCTION 

A number of recent experiments have studied the be- 
havior of a small number iV (1 ~ 100) of electrons which 
have been confined to two-dimensional quantum dots. 
When these quantum dots are subjected to large per- 
pendicular magnetic fields, a rich phenomenology may 
be observed, much of it due to the same physics under- 
lying the bulk quantum Hall effect 0J^. In particular, 
there has been a great deal of interest in the maximum 
density droplet (MDD) [|-|l6|, which is the small TV ana- 
log of the bulk v = 1 state. Experimentally, the MDD 
has been identified in several studies, including those of 
Ashoori et al. ^ , Klein et al. |^,^, and Oosterkamp et 
al 1^. Because the MDD's size decreases monotonically 
with magnetic field, it is clear from electrostatic consid- 
erations that a sufficiently strong magnetic field must 
destabilize the MDD into a lower density state. A num- 
ber of theoretical studies have discussed a type of insta- 
bility called an edge reconstruction, in which the decrease 
in charge density first occurs along a ring just within 
the MDD's perimeter. In this article we demonstrate 
that, neglecting the possibility of spin texturing |@-^, 
the MDD destabilizes into a state in which the droplet 
edge is best described as a one-dimensional Wigner crys- 
tal. The edge reconstructed density profile may then be 
understood as an angular momentum projection of the 
crystallized edge. 

In an early study of spin-polarized edge reconstruc- 
tions, Chamon and Wen |10|] identified approximate 
ground states of the quantum dot by minimizing its en- 
ergy over single Slater determinants of symmetric gauge 
lowest Landau level states. They found that increasing 
the magnetic field forces a transition in which several 
electrons simultaneously move from adjacent states in the 
interior of the droplet to states at the edge. This state, 
hereafter referred to as the compact ring (CR) state, 
therefore consists oi a, v = 1 central disk surrounded by 



a ring in which all single particle states are completely 
occupied. 

Other studies on smaller dots have revealed a rather 
different structure. For example, Maksym has stud- 
ied the electron-electron correlation functions of small 
droplets {N < 6) using exact diagonalization |l^]. His 
work finds that for v the electrons occupy distinct, 
localized positions relative to each other in a structure 
which resembles a small piece of a Wigner crystal [ p^ . 
In a study of somewhat larger dots, {N < 20), MuUer 
and Koonin [O] found that for ^ 1, an approximate 
Hartree-Fock (HF) Hamiltonian is minimized by a state 
whose density profile consists of a ^ 1 central droplet 
surrounded by a ring of localized electrons. 

Given that small dots are expected to Wigner crystal- 
lize at high magnetic fields, it seems reasonable to expect 
that similar behavior would occur in larger dots. How- 
ever, the crystallization of large dots would differ from 
that of small dots in at least one important respect: it 
would occur nonuniformly. This is most evident if one 
considers a droplet which begins as an MDD, and then 
evolves into a crystallized state as the magnetic field in- 
creases. In this case, any portion of the dot which crystal- 
lizes must be partially depopulated, and the extra charge 
moved to orbitals beyond the MDD perimeter. Since this 
transfer of charge costs less confinement energy for elec- 
trons which begin near the droplet edge, we expect the 
droplet edge to crystallize first, while the droplet center 
initially remains liquid | p9[ . 

In the following discussion, we will present numerical 
evidence that the MDD does, in fact, destabilize into an 
edge crystallized state when the magnetic field increases. 
To be precise, we will present a partial diagonalization 
of the many body Hamiltonian with the deep interior 
orbitals frozen. We will then calculate the overlap of 
this ground state with the "projected necklace" (PN), 
a trial wavefunction which exhibits edge crystallization. 
See Fig. (2). As a result, we find that PN states have 



1 



high overlap with diagonalized ground states, and have 
lower energy than Chamon and Wen's compact ring (CR) 
state. In addition, we have found that a reconstruction 
to a PN state may exhibit a much weaker addition spec- 
trum signature than does a reconstruction to a CR state. 
Whether the reconstruction signature is, in fact, weak 
depends on the detailed form of the confinement poten- 
tial. 



II. PROJECTED WIGNER NECKLACE STATES 

We begin our study by considering parabolically con- 
fined two-dimensional electrons in the presence of a per- 
pendicular magnetic field. The noninteracting Hamilto- 
nian for this problem is 



(1) 



where m* is the band mass, loq is the confinement 
strength (typically fiUo ~ 1-3 meV), and A = ^(z x r). 
For sufficiently high B, the only appreciably occupied 
eigenstates of Hq are the wavefunctions of the lowest 
Fock-Darwin level (LFDL) g, 

Vz(r) = (2'+Vn)-V2(l)'e-l-lV4/^, (2) 

where / is the angular momentum, z = x — iy, and Iq 
is the characteristic length, given by Iq — lmK~^, with 
= , K = (1 + 44)^/^ and ujc = The 

angular momentum of ijji in I, and its energy is 



ei = hiLOc {k^ + 1) - I) 



(3) 



Using the ipi^s as our basis states, we introduce the 
second-quantized Hamiltonian for many-electron quan- 
tum dots. 



(4) 



tjmn 



Here, Vijmn is the Coulomb interaction |g 



magnetic field, its general inward motion will be punc- 
tuated by a series of transitions in which the dot jumps 
to ground states of larger area, and accordingly, larger 
angular momentum, in order to decrease its Coulomb 
energy. 

The state of lowest total angular momentum L en- 
tirely contained in the spin-polarized LFDL is the MDD 
|]lO| , P^ , which is simply a single Slater determinant of 
the N LFDL single particle states with angular momenta 
0, 1, . . . , A^— 1. Consequently, when N is small, the MDD 
is the first state in the spin-polarized LFDL to become 
a ground state of the quantum dot as B increases from 
zero. A further increase in B will compress the MDD, so 
that it eventually becomes smaller in diameter than the 
classical electrostatic solution for a parabolically confined 
2D droplet, which is a hemisphere |^l|. At this point the 
contrast between the electronic density and the electro- 
static solution is most pronounced at the droplet edge: 
The electronic density slightly inside (outside) the MDD 
edge is greater (less) than would be expected from the 
electrostatic solution. Consequently, the initial destabi- 
lization of the MDD at high B involves a redistribution 
of the charge density at the MDD edge, and leaves the 
V = \ center intact |^ . 

We seek to characterize these instabilities by introduc- 
ing a set of trial wavefunctions which incorporate the 
possibility of crystalline correlations at the v — 1 edge. 
To begin, we construct a state which consists of a com- 
pact central disk of D electrons surrounded by a necklace 
oi R = N — D localized electrons equally spaced along 
a ring of radius u. Electrons in the central disk are as- 
sumed to occupy wavefunctions tpi, given by Eq. (H). 
Ring electrons are assumed to have orbitals 



r2 



and c|(q) is the creation (annihilation) operator corre- 
sponding to tpi- 

For the special case of parabolic confinement, the 
eigenstates of H depend on B only through the Iq depen- 
dence of the single particle orbitals. When B increases, Iq 
decreases, and the many-electron states squeeze inward, 
thereby increasing their Coulomb energy. It follows that, 
when a quantum dot is subjected to a steadily increasing 



Xp(r) = (27r)-i/2 



exp 



-^(|r-Upp-2zz-r xup) 



(6) 



where Up — u{cos4>p, sin(/>p) are the guiding centers, and 

R 



f>p — with p = 0, . . . , i? — 1. Our many-body wave- 



function is thus 
{rQ,...,rN-i\N,R,u;WN) [] M^k) H Xi~D{ri) 



'D-\ N-1 



, fc=0 l=D 



[V'm(ri)V'n(r2), 

(5) 



(7) 



where A is the antisymmetrization operator. We call 
this state a "Wigner necklace" . The density profile of a 
Wigner necklace with N = iO, R = 18, and u = 10.5/o is 
plotted in Fig. (1). 
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FIG. 1. Density contours of a Wigner necklace state with 
N = 40, 18, and u = lO.SZo- 

In order to project | A^, R, u; WN) onto fixed total an- 
gular momentum L, we seek to expand it into a sum over 
Slater determinants of the ipi^s. Because only states with 
a D electron compact central disk contribute to the ex- 
pansion, it is convenient to define |D, toi, . . . , mji) as an 
N body state consisting of a I? electron compact central 
disk and R additional single particle orbitals, that is, 



(ro, 



,rN-i I D,mi, 



,mR) 



(D-l N-l \ 

n M^k) n v^m,-x,+.(ro 
k=0 l=D ) 



(8) 



Defining P(V) to be the operator which projects onto 
total angular momentum L, we thus have 

P{L)\N,R,u-WN)= 

mx + ...+mR=L-\D(D-\) 

ID, mi, ... , mR){p, mi, ... , mi^|7V, i?, u; WN), (9) 

where the sum only includes Slater determinants with 
the correct total angular momentum. In practice, it 
is necessary to truncate this sum by disposing of all 
Slater determinants involving single particle momenta 
with rrii > N + k — 1 + W, where W acts as a cutoff. 
In this work, we use W >7. The coefficients in Eq. (|), 
which are computed in the Appendix, are given by 



{D,mi,...,mR\N,R,u;WN) = 

R , 

C^{W)-^I^ n sin^(m,-m,) 



\<i<j<R 



(10) 



where C is a constant common to every coefficient in the 
expansion of a particular |iV, R, u; WN). 



It follows from (see the Appendix) that if we define 
P(V) to be the operator which projects states onto total 
angular momentum L, then P(V)\N ,R,u\W N) = un- 
less L = - 1) -I- fci?, where fc = 0, 1, . . . (A similar 
result has been derived in Refs. and Q for the 
case of totally crystallized droplets.). Hence, the angular 
momentum projected states are conveniently defined as 

|iV, i?, fc; PN) = P{^N{N - 1) + kR)\N, R, u; WN), 

(11) 

where PN means "projected necklace." Note that 
\N, R, k; PN) is independent of the ring radius u used 
to construct the unprojected \N, R,u;WN) . The k — 1 
states are particle- hole excitations of the MDD . 

To illustrate the effect of the angular momentum 
projection, we examine the electron-electron correlation 
function 



P(r,ri) = (*t(r)^ft (1.^)^(1.^)^(1.))^ 



(12) 



where *(r) = T,ZoMr)ci Fig- (2)> we plot 

P{r, ri) for the state |40, 18, 3, PN), where we have fixed 
ri = {x,y) = (7.5/o,0). 
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FIG. 2. Two particle correlation function P(r, ri) for the 
state |40, 18, 3; PN), plotted as a function of r, with ri fixed 
at ri = (7.5^0,0). P(r, ri) was rescaled so that its maximum 
value is one. In order to emphasize the state's axial modula- 
tion, extra contours were added at levels 0.855 and 0.865. 

We note that this state is the lowest energy PN state 
with iV = 40 when hcu = 1.6meV and B = 2.42T. See 
Fig. (6) and the accompanying discussion, which will 
follow. Indeed, a modulation of P(r, ri) is evident at the 
droplet's edge, but is much weaker than the localization 
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in the unprojected WN. Extra contours were added to 
the figure in order to reveal the effect (see the caption). 

It is instructive to note the similarities between the PN 
states and the compact ring (CR) states discussed ear- 
lier. Formally, we define the CR state |iV, R, k; CR) to be 
a Slater determinant constructed from the single particle 
LFDL states with angular momenta I — 0, . . . , D — 1 and 
angular momenta I — D + k,...,N + k — 1, where D ~ 
N -R, as above. Both | A^, R, k; CR) and |iV, R, k; PN) 
may be thought of as states formed by moving k elec- 
trons from the bulk to the edge of an N electron MDD 
in order to separate an R electron ring from an N — R 
electron core. Furthermore, for a given triplet {N,R, k), 
\N, R, k; CR) and | A'', R, k; PN) have the same total an- 
gular momentum L. It follows from Eqs. (||) and (jj) that 
the sum of kinetic and confinement energies of the two 
states is the same for all B. Therefore, the state with the 
lower Coulomb energy will be favored for all B. 

Motivated by this observation, we have made pair- 
wise comparisons of the Coulomb energies of the states 
\N,R,k;CR) and \N,R,k;PN) for all {N,R,k) such 
that 25 < iV < 60, 2 < i? < 25, and 1 < fc < 10. The 
ranges of R and k have been selected to include all PN 
and CR states which are likely to become ground states 
of the quantum dot at magnetic fields slightly beyond 
the stability range of the MDD. We have confirmed this 
selection with appropriate minimizations of the droplet 
energy with respect to PN states, and in separate calcu- 
lations, with respect to CR states. 

For nearly each triplet {N, R, k) within the stated 
range, we find that | A^, R, k; PN) has the lower Coulomb 
energy, although the difference is small, typically around 
0.1% '^0.3%. A sample of this comparison is given in 
Table (1). 





Compact 


Projected 




{N,R,k) 


Ring 


Necklace 


E{CR) - E{PN) 


(40,18,1) 


124.577 


124.487 


0.090 


(40,18,2) 


123.350 


123.222 


0.128 


(40,18,3) 


122.123 


121.982 


0.141 


(40,18,4) 


120.914 


120.771 


0.143 


(40,18,5) 


119.732 


119.593 


0.139 


(40,18,6) 


118.585 


118.450 


0.134 


(40,18,7) 


117.475 


117.345 


0.131 


(40,18,8) 


116.406 


116.277 


0.129 


(40,18,9) 


115.377 


115.247 


0.130 


(40,18,10) 


114.389 


114.260 


0.129 



TABLE I. Comparison of the Coulomb energies of com- 
pact ring states and projected necklace states with A*' = 40, 

2 

R = 18, and fe = 3. All energies are in units of 4- . 



The only triplets in the stated range for which the 
corresponding CR state has the lower Coulomb energy 
are {N, 2, 10) for A^ = 25, . . . , 30 and (25, 3, 10). However, 
these exceptions are not a concern because in this range 
of A^, i? ^ 3 for the ground state, regardless of whether 
one assumes that the ground states are CR states or PN 
states. These results are encouraging, as they show that 
in the vicinity of the edge reconstruction, the ground 
state of a HF calculation will necessarily be unstable with 
respect to at least one PN state. 



III. DIAGONALIZATION CALCULATIONS 

We wish to further test the validity of the PN states 
by comparing them with the results of diagonalization 
calculations. Unfortunately, the systems we are studying 
are too large to accommodate complete diagonalization. 
Instead, we restrict the basis of our diagonalization cal- 
culation to Slater determinants of LFDL states which 
include a compact central disk of D electrons, defin- 
ing \N, R, k; CS) to be the lowest energy eigenstate of 
H within this basis for which L = Lmdd + kR, where 
Lmdd = ^N{N — 1). This restriction is reasonable be- 
cause, as discussed above, ground states with L > Lmdd 
are expected to retain a. u ^ 1 center. 

In Fig. (3), we display the single particle occu- 
pancies of the states |40, 18, 3, PA^), |40, 18, 3, Ci?), and 
|40, 18,3,GS'). 
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FIG. 3. Panels (a), (b), and (c) show, respectively, the 
single particle occupancies of the projected necklace state 
|40, 18, 3; PN), compact ring state |40, 18, 3; CR), and diago- 
nalized state |40, 18, 3; GS). 

Clearly, the PN state resembles the diagonalized state 
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more closely than the CR state does. The main differ- 
ence between the PN state and the diagonalized state is 
that the occupancy of the PN state drops sharply at the 
edge of the central disk. This discrepancy becomes more 
pronounced as the separation of the ring and central disk 
increases. 

More quan- 
titatively, we have computed (N, R,k; PN\N, R,k;GS), 
the overlap between the diagonalized ground state and 
the projected necklace state, for states with = 30 and 
N = 40. The results are plotted in Fig. (4). 
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FIG. 4. Overlaps of projected necklace states with ground 
states of diagonalization calculations. In panels (a) and (b), 
A'' = 30 and A*' = 40 electrons, respectively. The data in both 
panels has the pattern of descending staircases, with each step 
corresponding to a single value of k. 

Interestingly, {N, R, k; PN\N, R, k; GS) is determined 
primarily by k. This is evident in Fig. (4), where the 
data points assemble into a staircase of neat horizontal 
bunches which descend with increasing k. The numeri- 
cal values of the overlaps presented range between 0.8~ 
I.G with the higher overlaps corresponding to smaller 
k. The k — 1 states, whose overlaps consistently ex- 
ceed 0.99, are particularly successful [^. On the other 
hand, {N,R,k;PN\N,R,k;GS) often drops to nearly 
zero (i.e. < 0.002) with increasing k. This effect oc- 
curs when some of the ring electrons in the diagonalized 
ground state become adjoined to the central disk, so that 
|A^, R, fc; GS) has, in effect, fewer electrons in its outer 
ring than \N, R,k; PN). As A^ increases, the preferred 
number of electrons in the outer ring increases, and this 
effect disappears for fixed {R,k). These results, along 
with the energetic studies discussed above, demonstrate 
that the PN states with small k are good approximations 



to the edge reconstructed states. 

We have also computed (A^, R, k; GR\N, R, k; GS), the 
overlap of a compact ring state with the corresponding di- 
agonalized state, for the same triplets (A^, R, k) for which 
we computed (A^, R, k; PN\N, R, k; GS). The results are 
plotted in Fig. (5). 



06 

Hi" 



(a) 





R=12 





R=14 




H=16 


□ 


R=18 


B 


R=20 



(b) 









O R=14 






A H=16 




7 


□ R=18 






a R=20 










06 


□ 






A 


A 


05 




6_ 








04 


X 







002/^ 



A 



0.02 




40 



FIG. 5. Overlaps of compact ring (CR) states with ground 
states of diagonalization calculations. In panels (a) and (b) 
N — 30 and A'^ = 40 electrons, respectively. The CR states 
are less accurate for the larger N, and are consistently less 
accurate than the projected necklace states. 

The computed values of {N, R, k; GR\N, R, k; GS) fall 
in the range 0.3-0.8, except for a few values which 
are nearly zero (i.e. < 0.002). The set of triplets 
{N,R,k) for which {N, R, k; CR\N, R, k; GS) is nearly 
zero is exactly the same set of triplets for which 
{N, R, k; PN\N, R, k; GS) is nearly zero. As for the ma- 
jority of triplets, for which the overlap is substantial 
(> 0.3), (N,R,k;GR\N,R,k]GS) is on average 36% 
less than the corresponding PN overlap when A^ = 30, 
and 44% less than the corresponding PN overlap when 
A^ = 40. 



IV. ADDITION SPECTRUM CALCULATIONS 

We now discuss the consequences of the PN states for 
an important experimental quantity, the quantum dot 
addition spectrum. An addition spectrum consists of a 
sequence of measurements of the chemical potential 



^,{N)=EQiN)-EoiN-l), 



(13) 



over some range of A^, where Eo{N) is the ground state 
energy of an A^-electron quantum dot M-M . To perform 
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this computation, we assume that the ground state in the 
1/ ^ 1 regime is the PN state which minimizes the energy, 
that is, we take 

EoiN) =^mm{N,R,k;PN\H\N,R,k;PN). (14) 

R,k 

The B dependence of /i(iV) for a parabohcally confined 
droplet with huj = 1.6meV, as computed from (13) and 
(|l|), is plotted in Fig. (6) for = 40 to iV = 46. 




FIG. 6. Addition spectrum traces for a droplet subject to a 
parabolic confinement of strength Hljo = l.GmeV, computed 
by energy minimization over projected necklace states. The 
A'^th line (as marked) corresponds to the addition of the TVth 
electron into the droplet. Each arrow indicates a unit increase 
in the A*' electron droplet's ground state value of k. The left- 
most column of arrows denotes the edge reconstruction. 



of Ref. |lO| , in which the onset of the reconstruction in- 
volves the simultaneous transfer of several electrons, and 
causes a sharp upward jump in the addition spectrum. 
However, because our approach is based on a more ac- 
curate solution of the same Hamiltonian (see prece ding 
discussion), we believe that the cusps seen in Ref. |lO|] 
are an artifact of the HF approximation. 

On the other hand, experimental work has demon- 
strated that the high B destabilization of the MDD 
causes distinct cusps in the addition spectrum of the 
quantum dot [^J^. In light of the results presented in 
Fig. (6), it is thus reasonable to ask whether circum- 
stances exist in which the experimental observation of 
addition spectrum cusps would even be consistent with 
spin polarized edge reconstructions. 

In fact, we may reintroduce edge reconstruction in- 
duced cusps into the addition spectrum by introducing a 
non-parabolic term into the confinement potential. To be 
specific, we write V{r) as the sum of the parabolic con- 
finement potential and an additional coffee cup shaped 
potential, giving 

V{r) = ^m*Lu'^r'^ + (3{r-a)9{r-a), (15) 

where /? and a are constants, and 9{x) is the Heaviside 
step function. The second quantized many-body Hamil- 
tonian now becomes 



E 

ijmn 



ijmn 



(16) 



with ji = (3 /^°° cPr^j*{v){r - 0)1/1; (rV 

By exact diagonalization of Eq. (16) within the LFDL, 
we compute the addition spectra of dots subject to and 
not subject to a coffee-cup confinement, and plot the re- 
sults, respectively, in Figs. (7) and (8). 



Notably, the addition spectrum traces are nearly fea- 
tureless in the vicinity of the edge reconstruction, except 
for small kinks due to increments in k. Recall that k 
may be understood as the number of holes separating 
the central disk and the compact edge. As B increases, 
the ground state value of k for a given N always increases 
monotonically in unit increments. The first increment of 
k for a given iV, from fc = to fc = 1, initiates the edge 
reconstruction. Each increment in the ground state value 
of k in the N electron dot is indicated in Fig. (6) by an 
arrow pointing to the A'^th addition line. The ground 
state value of R at the onset of the edge reconstruction 
is determined mainly by the requirement that the fill- 
ing fraction at the center of the reconstructed edge be 
V 1. As N increases from A^ = 30 to A^ = 49, this 
value increases from i? = 16 to i? = 20. 

These results are in distinct contrast with the HF result 
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FIG. 7. Addition spectrum traces, computed with exact di- 
agonalization, for a droplet subject to both a parabolic con- 
finement potential of strength huio = S.OmeV, and a non- 
parabolic confinement potential characterized by parameters 
a = 40nm, and /3 = 2.0meV/nm (See Eq. (Ilq)). Compare 
with Fig. H. 



(7), a = 40nm and f3 = 2.0meV/nm. In these plots, a 
lone arrow near each addition spectrum trace marks the 
edge reconstruction. Indeed, the edge reconstruction in- 
duced kinks of the nonparabolically confined system are 
much larger than those of the parabolically confined sys- 
tem. The kink size also increases with electron number, 
similar to the data of Oosterkamp et al ||] . 

In summary, we have shown that projected necklace 
(PN) states accurately describe the spin-polarized insta- 
bilities of 1/ 1 quantum Hall droplets. PN states whose 
total angular momentum differs from that of the MDD 
by 70h or less typically have overlaps with diagonalized 
states in the range 0.8-0.99. By comparison, the com- 
pact ring (CR) states, which are a generalization of the 
states identified by Chamon and Wen [0 , have overlaps 
of 0.3-0.8 with diagonalized states. We have also shown 
that near the edge reconstruction, the lowest energy PN 
state is always lower in energy than the lowest energy 
CR state. By performing energy minimizations over PN 
states, we have shown that edge reconstruction induced 
cusps in the addition spectra of parabolically confined 
dots are much smaller than suggested by Hartree-Fock 
JlOf studies. Finally, we have shown that large cusps in 
the addition spectra may occur in dots with nonparabolic 
confinement potentials. 
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FIG. 8. Addition spectrum traces, computed with exact 
diagonalization, for a droplet subject to a parabolic confine- 
ment of strength hcoo — S.OmeV. Compare with Fig. |^ 

For both dots, hujQ = S.OmeV, and for the dot of Fig. 
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APPENDIX: 

In this appendix, we first compute the coefficient 
{D,mi, . . . ,mii\N,R,u;WN), introduced in Eq. (|). 
Recall that in general, the overlap of two Slater deter- 
minants |ai, . . . , qk) and |5i, . . . , bx), where the Ui and 
bi are single particle states, is 



[ai, . . 

Hence 

{D,mi, 



,aK\bi 



, bjA = 



(ai|6i) 
{ailbx) 



(flA'l^l) 

{aKlbx) 



,mR\N,R,u;WN) = 



Id Odxr 
Ob.xd Sfl 



(Al) 



(A2) 



where Id is the D x D unit matrix, due to the inclusion 



of the D electron compact disk in each state, 0^ 
m X n zero matrix, and 



Sr = 




{fprnnlXR-l) 



is the 



(A3) 
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To compute S/j, we use the result (V'ilXp) = aze*''^", 
where ai = (2'?!)~-^/^M'e~" Z'*. Factoring the q;;'s out 
of the determinant leaves in the form of an alternant, 
which may be readily evaluated, yielding 

mi, . . . , TOfllTV, R, u- WN) = 

ft 

CYlim^l)-'^' n sin(^(m,-m,)Y (A4) 

1=1 l<i<j<R ^ ^ 

where C is a constant common to each term in the expan- 
sion of a particular |iV, R, u; WN). Because of the factor 

Yli^j sm(^jj {nii ~mj)^, {D,mi, . . . ,mR\N, R,u;WN) 

must vanish if mi = rrij mod R for any pair {mi^rrij). 
This result, along with the observation that there are 
exactly R m-i's, lead us to the following lemma: 

Lemma 1 For each integer j,D<j<D + R — 1, each 
state \D, mi, . . . , mji) in the expansion of \N, R, u; WN) 
must include exactly one single particle state (ex- 
cluding states in the compact central disk) for which 
/ij = J mod i?. Specifically, we write /ij — j = kjR, 
where kj is a non- negative integer. 

It follows that the total angular momentum of the 
states in the ring differs from the minimum possible value 
^R{2D + R~ 1) hy kR, where k = J^f^j^^^ kj. The 
significance of k is discussed in the main body of the 
text. The lemma is very useful for computation, as it 
greatly reduces the number of terms in the expansion 
of a state |iV, i?, fc; PN), and allows us to project wave- 
functions \N,R^u\WN) onto values of L which would 
otherwise be inaccessible. 
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